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We present the first direct measurement of three-body interactions in a colloidal system comprised 
of three charged colloidal particles. Two of the particles have been confined by means of a scanned 
laser tweezers to a line-shaped optical trap where they diffused due to thermal fiuctuations. Upon 
the approach of a third particle, attractive three-body interactions have been observed. The results 
are in qualitative agreement with additionally performed nonlinear Poissson-Boltzmann calculations, 
which also allow us to investigate the microionic density distributions in the neighborhood of the 
interacting colloidal particles. 

a 

I. INTRODUCTION 

vn 

Pair interactions in dense systems are in general affected by the presence of many other surrounding particles. To 
,_i take such many-body interactions into account, the degrees of freedom of other particles are often integrated out, 
leading to effective pair potentials. This concept is often the only way to handle systems where a large number of 
O different length and time scales coexist. It is important to realize, however, that effective potentials - in contrast to 
^ true pair potentials - can not be regarded as fundamental quantities because their parameters depend on the state 
of the system. In addition, no unique way to derive the effective potentials exists and the effective pair potential 
picture very often leads to thermodynamic inconsistencies [T. Accordingly, a correct description of any liquid or 
solid must explicitly take into account many-body effects (and in particular three-body effects as the leading term). 
Already in 1943 it has been supposed by Axilrod and Teller (AT) [2] and later also by Barker and Henderson [3] 
(~| that three-body interactions may significantly contribute to the total interaction energy in noble gas systems. This 
O seems to be surprising because noble gas atoms posses a closed-shell electronic structure and are therefore often (and 
, ^ , erroneously) regarded as an example of a simple liquid . The conjecture of Axilrod and Teller, however, was confirmed 
only very recently, when large-scale molecular dynamics simulations for liquid xenon and krypton [U [S] was compared 
T-H with structure factor measurements at small q- vectors performed with small-angle neutron scattering [HlEj- In these 
^ papers it has clearly been demonstrated that only a combination of pair-potentials and three-body interactions, the 
latter in the form of the AT-triple-dipole term leads to a satisfactory agreement with the experimental data. In 
the meantime, it has been realized that many-body interactions have to be considered also for nuclear interactions 
[F, inter atomic potentials, electron screening in metals [5^, photo- ionization, island distribution on surfaces [10], and 
even for the simplest chemical processes in solids [12^ like breaking or making of a bond. 

In view of the general importance of many-body effects it seems surprising that until now no direct measurements of 
these interactions have been performed. This is largely due to the fact that in atomic systems, positional information 
is typically provided by structure factors or pair-correlation functions, i.e. in an integrated form. Direct measurements 
• • of many-body interactions, however, require direct positional information beyond the level of pair-correlations, which 
is not accessible in atomic or nuclear systems. In contrast to that, owing to the convenient time and length scales 
^> involved, the microscopic information is directly accessible in colloidal suspensions. In addition, the pair-interactions 
^ in colloidal suspensions can be varied over large ranges, e.g. from short-ranged steric to long-ranged electrostatic or 
even dipole-dipole interactions. 

In the present study we used charged colloidal particles whose interactions are mediated by the microscopic ions in 
the electrolyte. The pair interaction in such systems is directly related to the overlap of the ion clouds (double-layers) 
which form around the individual colloids and whose thickness is determined by the ionic strength of the solution. 
In highly de-ionized solutions, these double-layers can extend over considerable distances. If more than two colloids 
are close enough to be within the range of such an extended double layer, many-body interactions are inevitably the 
consequence. Accordingly, deviations from pair-wise additive interaction energies are expected in charge-stabilized 
colloidal systems under low salt conditions. 
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Here we present the first direct measurement of three-body interactions, performed in a suspension of charged 
colloidal particles. This was achieved by scanned optical tweezers, which provided a trapping potential for two colloidal 
particles. When a third particle was present, considerable deviations from pair-wise additive particle interactions have 
been observed. These deviations increased, as the distance of the third particle was decreased and were used to extract 
three-body interaction potentials. We have additionally performed non-linear Poisson-Boltzmann calculations for the 
same parameters and same configurations as chosen in the experiment. Deriving the interaction potentials from the 
solutions of the Poisson-Boltzmann equation, we have correctly taken three-body terms into account. The numerically 
obtained three-body potentials are in qualitative agreement with the experimental results. 

Experimental evidence for many-body interactions has been already obtained from effective pair-interaction poten- 
tial measurements of two-dimensional (2D) colloidal systems. Upon variation of the particle density, a characteristic 
dependence of the effective pair interaction was found which has been interpreted in terms of many-body interactions 
[13j . However, during those studies the relative contributions of different many-body terms could not be further 
resolved. Performing the experiment described in this paper, i.e. observing the system of only three particles, we 
were able to measure the three-body interactions directly. 

II. EXPERIMENTAL SYSTEM 

As colloidal particles we used charge-stabilized silica spheres with 990nm diameter suspended in water. A highly 
diluted suspension was confined in a silica glass cuvette with 200/im spacing. The cuvette was connected to a closed 
circuit, to deionize the suspension and thus to increase the interaction range between the spheres. This circuit 
consisted of the sample cell, an electrical conductivity meter, a vessel of ion exchange resin, a reservoir basin and a 
peristaltic pump |T3]. Before each measurement the water was pumped through the ion exchanger and typical ionic 
conductivities below 0.07/zS/cm were obtained. Afterwards a highly diluted colloidal suspension was injected into 
the cell, which was then disconnected from the circuit during the measurements. This procedure yielded stable and 
reproducible ionic conditions during the experiments. Due to the ion diffusion into the sample cell, the screening 
length decreased linearly with time during the measurements. The rate of change of the screening length, 
however, was only less than half a percent per hour, which means that in the time needed to perform a complete set 
of measurements, the ionic concentration did not change more than about 1 percent. This tiny variation has been 
taken into account when performing the PB calculations (see section IV.). 

First, three particles were brought in the field of view of the microscope after they had sedimented down to the 
bottom plate of the sample cell (FigjT|). Two particles were trapped with line-scanned optical tweezers, which was 
created by the beam of an argon ion laser being defiected by a computer-controlled galvanostatically driven mirror 
with a frequency of approximately 350Hz. The time averaged intensity along the scanned line was chosen to be 
Gaussian-distributed with the half-width ~ 4.5 fim. The laser intensity distribution perpendicular to the trap 
was given by the spot size of the laser focus, which is also Gaussian with ay ~ 0.5 /xm. This yielded an external 
laser potential acting as a stable quasistatic trap for the particles. Due to the negatively charged silica substrate, the 
particles also experience a repulsive vertical force, which is balanced by the particle weight and the vertical component 
of the light force. The potential in the vertical direction is much steeper than the in-plane laser potential, therefore 
vertical particle fiuctuations can be disregarded. The particles were imaged with a long-distance, high numerical 
aperture microscope objective (magnification x63) onto a CCD camera and the images were stored every 120 ms. 
The lateral positions of the particle centers were determined with a resolution of about 25 nm by a particle recognition 
algorithm. 

Three-body interaction potentials were measured in this setup by performing the following steps (which will be 
explained in detail below): First only one particle was inserted into the trap and its position probability distribution 
was evaluated from the recorded positions. From this the external laser potential ul could be extracted. Next, we 
inserted two particles in the trap and measured their distance distribution. From this, the pair-interaction potential 
was obtained. Finally, a third particle was made to approach to the optical trap by means of additional point optical 
tweezers (focus size ~ 1.3 /im), which held this particle at a fixed position during the measurement. From the distance 
distribution of the first two particles we obtained the total interaction potential for the three particles. Finally, we 
substracted a superposition of pair potentials (known from the previous two-particle measurements) from the total 
interaction energy to obtain the three-body interaction. 
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III. DATA EVALUATION AND EXPERIMENTAL RESULTS 



We first determined the external potential acting on a single particle due to the optical line trap. The probability 
distribution P{x, y) of finding a particle at the position (x, y) in the trap was evaluated from the recorded positions. 
P{x,y) depends only on the temperature and the external potential UL{x,y) created by the laser tweezers. According 
to the Boltzmann probability distribution P{x, y) = Pj^e~'^^^^-^'y\ with being a normalization constant. Taking the 
logarithm of P{x, y) yields the external potential ul{x, y) with an offset given by log P^. The probability distributions 
in X and y directions are statistically independent, and can therefore be factorized. The laser potential is thus 
ul{x, y) — UL{x) + UL{y). The potential along the x axis is shown in Figj2]for various laser intensities. As can be seen, 
all renormalized potentials fall, within our experimental resolution, on top of each other. This clearly demonstrates 
that the optical forces exerted on the particles scale linearly with the input laser intensity. This fact allows us to use 
different external laser powers for two-body and three-body experiments (in the three-body experiment, due to the 
additional repulsion of the third particle, a stronger laser power is needed to keep the mean distance between the two 
particles similar). The corresponding potential in the perpendicular {y) direction has the same (Gaussian) shape, but 
it is much steeper due to the chosen scanning direction. Therefore, the particles hardly move in the y direction during 
a measurement. 

Next, we inserted a second particle in the trap. The four-dimensional probability distribution is now 
P{xi,yi,X2,y2) = Pi2e~^^^^'--^^''^^'>^^^''^^'-^^^~^^'^'-''^\ with Xi and yi being the positions of the i-th particle relative 
to the laser potential minimum and U{r) the distance dependent pair-interaction potential between the particles. 
This can be projected to 



P{r) = J J J J P{xi,yi,X2,y2)5{yJ {xi - X2Y + (2/1 - 2/2)^ - dxidx2dyidy2 = 

= Piae-^U^-^) /"/"/"/" e-''("-(^i'y^)+"-(''='y^» ,5(v/(xi - xz)^ + (yi - y2)2 ~ r) dxidxzdyidys (1) 



In principle the integral is constituted of all possible configurations of two particles with distance r. Performing the full 
four-dimensional integration, however, is difficult because of the limited experimental statistics. This problem can be 
overcome by the following two considerations. First, due to the Gaussian shape of the external potential, the most likely 
particle configurations are symmetric with respect to the potential minimum of (any asymmetric configuration for 
constant r has a higher energy) . Secondly, particle displacements in ?/-direction are energetically unfavorable because 
ax 3> (Ty. Accordingly, for r ~ const the minimum energy configuration is {xi — r/2,yi — 0,X2 — — r/2,?/2 = 0). It 
has been confirmed by a simple calculation with the experimental parameters that all other configurations account 
for only for less than 1 percent of the value of the integral in Eq.([T]). Accordingly, Eq.Q reduces to 

P(r) = Ppe-''(u('^)+2"^('-/2.o)) . (2) 

Since ul{x, y) is known from the previous one-colloid measurement, we can obtain the interaction potential U{r) from 
the measured P{r), 

pU{r) = - log P(r) - 2/3ML(r/2,0) + logPo . (3) 

The normalization constant Pq was chosen in a way that U{r) — > for large particle separations r. We first measured 
U{r) according to the above procedure in the absence of a third particle. As expected, the negatively charged colloids 
experience a strong electrostatic repulsion which increases with decreasing distance. The pair-interaction potential of 
two charged spherical particles in the bulk is well known to be described by a Yukawa potential [151 ES] 



where Z* is the renormalized charge |17j of the particles, Xb the Bjerrum- length characterizing the solvent (As = 
/AneeokBT, with e the dielectric constant of the solvent and e the elementary charge), the Debye screening length 
(given by the salt concentration in the solution), R the particle radius and r the centre-centre distance of the particles. 
Fig j3] shows the experimentally determined pair-potential (symbols) together with a fit to Eq.Q (solid line). As can 
be seen, our data are well described by Eq.(|4]). As fitting parameters we obtained Z* « 6500 electron charges and 
Ki 470nm, respectively. The renormalized charge is in good agreement with the predicted value of the saturated 
effective charge of our particles [TH] [12] and the screening length agrees reasonably with the bulk salt concentration 
in our suspension as obtained from the ionic conductivity. Given the additional presence of a charged substrate, it 



4 



might seem surprising that Eq.Q describes our data successfully. However, it has been demonstrated experimentally 
[20] and theoretically [211 ES] that a Yukawa-potential captures the leading order interaction also for colloids close 
to a charged wall. A confining wall introduces only a very weak (below 0.1 fc^T) correction due to additional dipole 
repulsion. This correction is below our experimental resolution. Repeating the two-body measurements with different 
laser intensities (50mW to 600mW) yielded within our experimental resolution identical pair potential parameters. 
This also demonstrates that possible light-induced particle interactions (e.g. optical binding [23]) are neglegible. The 
approach of the third particle by means of an additional optical trap could, in principle, lead to additional light- 
induced interactions between the laser spot and the two particles kept in the line trap. To exclude such effects, 
we repeated the two-particle measurements and approached an empty trap (without the third particle) to the line 
trap where the two particles were fluctuating. Within our experimental resolution, we again observed identical pair 
potentials, which suggests, that the additional optical trap has no influence on the two particles in the line trap. 
When a third particle is present at a distance d along the perpendicular bisector of the scanned laser line (cf. inset 
of Fig ij, the total interaction energy U{r, d) is not simply given by the sum of the pair-interaction potentials Eq.Q 
alone but also contains an additional term. Following the definition of McMillan and Mayer [24] , U (r, d) is given by 



U{r,d) = Upair{ri2) + Upair{ri3) + Upair{r23) + Ui23{ri2,ri3,r23) , (5) 

with Upair{fij) being the pair-potential between particles i and j as defined in Eq.(j4]) and U123 the three-body 
interaction potential. Distances ri2, ^23 and ri3 are the distances between the three particles which can, due to 
the chosen symmetric configuration (r23 = ris), be expressed by the two variables r = ri2 and d = — (r/2)^. 

We have followed the same procedure as described above for the case of two particles. First, we have measured the 
probability distribution P(r; d) of the two particles in the laser trap with the third particle fixed at distance d from the 
trap. Taking the logarithm of P{r;d) we extracted the total interaction energy U{r,d) [25^. The results are plotted 
as symbols in Fig|4] for the distance of the third particle d — 4.1, 3.1, 2.5 and 1.6 /im, respectively. As expected, 
U{r, d) becomes larger as d decreases due to the additional repulsion between the two particles in the trap and the 
third particle. In order to test whether the interaction potential can be understood in terms of a pure superposition 
of pair-interactions, we first calculated U{r,d) according to Eq.([5]) with U123 = 0. This was easily achieved because 
the positions of all three particles were determined during the experiment and the distance-dependent pair-potential 
is known from the two-particle measurement described above (Figj3]). The results are plotted as dashed lines in Figj4] 
Considerable deviations from the experimental data can be observed, in particular at smaller d. These deviations can 
only be explained, if we take three-body interactions into account. Obviously, at the largest distance, i.e. d — 4.1/im 
our data are well described by a sum over pair-potentials which is not surprising, since the third particle cannot 
infiuence the interaction between the other two, if it is far away from both. In agreement with theoretical predictions 
[26], the three-body interactions therefore decrease with increasing distance d. 

According to Eq.(|5| the three-body interaction potential is simply given by the difference between the measured 
U{r, d) and the sum of the pair-potentials (i.e. by the difference between the measured data and their corresponding 
lines in Fig|4|. The results are plotted as symbols in Figjs] It is clearly seen, that in the case of charged colloids U123 
is entirely attractive and becomes stronger as the third particle approaches. It is also interesting to see that the range 
of M123 is of the same order as the pair-interaction potentials. It might seem surprising that it is possible to sample 
the potential up to energies of IS/csT, as conflgurations of such a high energy statistically happen only with very low 
probability. In this experiment we can choose the energetic range of the potential we want to sample by adjusting the 
strength of the line tweezers. The laser potential pushes the particles together, which allows us to sample different 
ranges of the electrostatic potential. Thus, to achieve a better resolution for smaller particle separations (e.g. higher 
potential values), the strength of the line tweezers had to be increased. The shape of the external potential was 
independent of the strength of the laser beam (see Figj2| and the magnitude scaled linearly with the input laser 
power. This allowed us to adjust the input laser intensity so as to obtain a suitable particle separation range. The 
external potential was obtained simply by scaling the Gaussian shown in Fig|2] 



IV. NUMERICAL CALCULATIONS 



In order to get more information about three-body potentials in colloidal systems, we additionally performed 
non-linear Poisson-Boltzmann (PB) calculations, in a similar way as in ^26j. The PB theory provides a mean-fleld 
description in which the micro-ions in the solvent are treated within a continuum approach, neglecting correlation 
effects between the micro-ions. It has repeatedly been demonstrated [37] US] that in case of monovalent micro- ions the 
PB theory provides a reliable description of colloidal interactions. The interactions among colloids arc, on this level, 
mediated by the continuous distribution of the microions and can be obtained once the local electrostatic potential 
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due to the microionic distribution is known. The normaHzed electrostatic potential ip{x,y,z), which is the solution 
of the non-linear PB equation, 

^'^'tjj^f) — K^sinh'0(r^ , 

n-'Vtp ~ AttXbu , r on colloid surface , (6) 

describes the equilibrium distribution of the microions around a given macroionic configuration. Here k is the inverse 
Debye screening length , the Bjerrum length (A^ = 0.72nm for aqueous solutions at room temperature) and a 
is the surface charge density on the colloid surface (constant charge boundaries are assumed for all colloids in the 
system), n is the normal unit vector on the colloid surface. We used the multi-centered technique, described and 
tested in other studies [29l [30] to solve the PB equation ^ at fixed configurations of three colloids and obtain the 
electrostatic potential tjj{x,y,z), which is related to the micro ionic charge density. Integrating the stress tensor, 
depending on ip{x,y,z), over a surface enclosing one particle, results in the force acting on this particle. First, we 
calculated how the force /12, and from it the pair-potential between two particles, depend on the distance between 
isolated two particles. Choosing the suitable bare charge on the colloid surface, we were able to reproduce the 
measured pair-interaction in FigjSj The calculation of three-body potentials was then carried out by calculating the 
total force acting on one particle in the line trap (say, particle 1) in the presence of all three particles and subtracting 
the corresponding pair- forces /12 and /13 obtained previously in the two-particle calculation. If there is any difference 
between the force on particle 1 obtained from the full PB solution for the three particle configuration and the sum 
of two two-body forces, this difference is due to the three-body interactions in the system. The difference is then 
integrated to obtain the three-body potential. The results are plotted as dashed and dotted lines in FigjS] and show 
qualitative agreement with the experimental data. To account for the deviations from the experimental data one has 
to take into account the following points: (i) there is a limited experimental accuracy to which the light potential can 
be determined. The accuracy decreases with increasing laser intensity (note that normalized potentials are plotted 
in Figj2]). In the three-body experiments, due to the presence of the third repulsive particle, a stronger light field is 
needed and the experimental error in determining the light potential is estimated to be around zblfcsT. Since we have 
to subtract the light potential twice from the total potential to obtain the three-body potential, this error doubles 
and we expect an error of about ±2kBT in the final result, (ii) An error of about -izlksT should be expected in 
the numerically obtained three-body potentials as well, (iii) While in the numerical calculation we assume identical 
colloidal spheres, in the experiment small differences with respect to the size and the surface charge are unavoidable. 
This effect, however, is rather small and leads to deviations on the order of 5 percent of the total potential, (iv) 
The numerical calculations do not take into account any effects which may be caused by the substrate. Although 
we expect such effects to be rather small (similar to the contribution for the pair potential) they can not completely 
ruled out. Considering the above mentioned uncertainties it should be emphasized that in particular the sign and the 
order of magnitude of the calculated potential compares well with our measured results. This strongly supports our 
interpretation of the experimental results in terms of the three-body interactions. 

We have measured and calculated the three-body interaction on a mesoscopic level, but since the colloidal interac- 
tions are mediated by the microions distributed in an electrolyte around the colloids, it is interesting to explore what 
happens on a microscopic level, i.e., what feature of the microscopic distributions lead to the observed three-body 
interactions. Of course, it is not possible to observe the microionic density experimentally, but in a Poisson-Boltzmann 
simulation such an information is easily accessible. Since the microion density depends monotonically on the elec- 
trostatic potential 'ip{r), it is enough to compare the electrostatic potentials to qualitatively discuss the microscopic 
picture. Of course, to a large extent, the potential ^{r) around three particles is just the superposition of potentials 
around individual particles, but since the solutions of nonlinear equations are in principle not superposable, we ex- 
pect to find small differences. It is indeed these small differences that are ultimately responsible for the three-body 
interaction. 

We started by reconsidering the two-particle problem. First, we solved the PB equation around a single isolated 
colloid to obtain the one-body non-linear electrostatic potential il]'^{r). Next we calculated the electrostatic potential 
ijp'if) for two colloidal particles at distance r and compared this potential to the superposition of two one-body 
potentials V'K'^ + "02 The difference is shown as a contour-plot in Fig|6^. It can be seen that micro-ions are 
rearranged in a complex way between the colloids. There is a weak additional polarization of the counter-ion cloud 
very close to the particle surfaces not captured by superposing the one-body potentials. However, all these effects are 
rather small and therefore, except for very small particle separations r, the superposed solution should still describe 
the two-body interactions with good accuracy. Not so for three particles. We have compared a superposition of three 
two-body electrostatic potentials with the correct non-linear three-body electrostatic potential [3T]. The difference is 
shown in Figj6]D. Obviously, differences are now much larger than in Figj6^. We notice that the counter-ion cloud 
polarization close to the colloid surface is correctly taken into account by two-body terms, while the ion distribution 
in the region among the colloids is poorly described by adding up two-body electrostatic potentials. There are fewer 
counter-ions in the region among the colloids than a pair- wise description predicts. This suggests that the entropy 
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gained by removing some excessive counterions (predicted by the superposition) from the inter-particle space is 
larger than the positive energy difference due to less efficient screening resulting from it. By integrating the potential 
difference from FigjBJ), one recovers the attractive three-body potential, already discussed, which is thus demonstrated 
to be a consequence of the nonlinearity of the physical equations governing the interactions in our system. The exact 
microscopic explanation of the phenomenon, however, is still lacking and further work is necessary to achieve it. 

V. CONCLUSIONS 

We have demonstrated that in case of three colloidal particles, three-body interactions are attractive and of the 
same range as pair-interactions. They present a considerable contribution to the total interaction energy and must 
inevitably be taken into account. Whenever dealing with systems comprised of many (i.e. more than three) particles, 
in principle also higher-order terms have to be considered. The relative weight of such higher-order terms depends 
on the particle number density p . While at low enough p a pure pair- wise description should be sufficient, with 
increasing density first three-body interactions and then higher-order terms come into play. We expect that there is 
an intermediate density regime, where the macroscopic properties of systems can be successfully described by taking 
into account only two- and three-body interactions [32]. Indeed liquid rare gases [S] and the island distribution of 
adsorbates on crystalline surfaces [10^ are examples where the thermodynamic properties are correctly captured by 
a description limited to pair- and thrcc-body interactions [33 . In colloidal systems we have shown the three-body 
interactions to be comparable in magnitude to the corresponding pair-interactions, therefore we there expect large 
macroscopic three-body effects in this intermediate density range. At even larger particle densities n-body terms with 
n > 3 have to be additionally considered, which may partially compensate. Even in this regime, however, many-body 
effects are not cancelled, but lead to notable effects, e.g. to a shift of the melting line in colloidal suspensions, as 
recently demonstrated by PB calculations [231 EHI • 

With some effort it is in principle possible to proceed to measure the higher order many-body terms and to study 
how the many-body expansion converges. Work on four-body interactions is in progress. 
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VII. FIGURE CAPTIONS 

Fig.l: Photograph of sample cell with two silica particles confined to a light trap created by an optical tweezers 
and a third particle trapped in a focused laser beam. The inset shows a schematic drawing of the experimental 
geometry. 

Fig. 2: The shape of the laser potential along the tweezers line for three different laser intensities (symbols: 
triangles lOOmW, circles 200mW and squares 500mW); for better comparison all curves are normalized to an 
intensity of lOOmW. The Gaussian fit is plotted as a solid line. 

Fig. 3: Measured pair-interaction potentials U{r) = Upair{r) (symbols) in the absence of the third particle. The 
data agree well with a DLVO potential, Eq.Q (solid line). In the inset the potential is multiplied by r and plotted 
logarithmically, so that the DLVO expression, Eq.Q transforms into a straight line. From a fit we obtained the 
effective charge Z* « 6500 and the screening length w 470nm. 

Fig. 4: Experimentally determined interaction energy U{r) (symbols) for two particles in a line tweezers in the 
presence of a fixed third particle with distance d on the perpendicular bisector of the line trap. For comparison the 
superposition of three pair-potentials is plotted as lines. Symbols and lines are labelled by the value of d. 

Fig. 5: Three-body potentials for different d. Measured three-body potentials indicated by symbols. The lines 
are three-body potentials as obtained from the solutions of the nonlinear Poisson-Boltzmann equation for three 
colloids arranged as in the experiment. The parameters in the Poisson-Boltzmann calculation were chosen so that 
the pair-interaction potentials were correctly reproduced. Symbols and lines are labeled by the value of d. 

Fig. 6: Contour plots of electrostatic potentials, a) Difference between the full electrostatic potential for two 
particles and the superposition of two one-particle potentials. The distance between the particles is r = 2.5/im. 
b) Difference between the full electrostatic potential for three particles and the superposition of three two-particle 
potentials. The distance between particle 1 and 2 is r = 2.5/xm and the position of the third particle is given by 
d = 1.6/im, being the closest distance realized in the experiments. The colour scales are in units of ksT. 
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